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Abstract 

We show that the one-dimensional extended Hubbard model has saturated ferro¬ 
magnetic ground states with the spin-triplet electron pair condensation in a certain 
range of parameters. The ground state wave functions with fixed electron numbers 
are explicitly obtained. We also construct two ground states in which both the spin- 
rotation and the gauge symmetries are broken, and show that these states are trans¬ 
ferred from one to the other by applying the edge operators. The edge operators are 
reduced to the Majorana fermions in a special case. These symmetry breaking ground 
states are shown to be stabilized by a superconducting mean field Hamiltonian which 
is related to the Kitaev chain with the charge-charge interaction. 
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1 Introduction 


The extended Hubbard model has been studied extensively to understand phenomena such as 
charge density wave, spin density wave and unconventional superconductivity which can not 
be described by the Hubbard model consisting of the electron hopping term and the on-site 
interaction term [DEIII3]. The Hamiltonian of the model is obtained by adding interaction 
terms of electrons on different sites to the Hubbard Hamiltonian. In the case where the 
added interaction together with the on-site one is dominant and is known to induce a certain 
ordering state with an energy gap, the model is well understood by considering the electron 
hopping as a perturbation. On the other hand, in order to understand phenomena which 
do not arise directly from interactions, we have to face the difficult problem of analyzing 
the interplay between the electron hopping and some interactions in a convincing way. The 
unconventional superconductivity corresponds to such a case. 

Here we restrict ourselves to the one-dimensional extended Hubbard model with nearest 
neighbour interactions. Despite the difficulty in analyzing correlated electron systems, there 
are a few rigorous results associated with superconductivity in this case. Most of these results 
are obtained through the Bethe ansatz method, and the superconducting ground states so 
far obtained are related to spin-singlet electron pair condensation 015116]. In this paper we 
provide another rigorous result for the model. By using a similar method in Ref. [7], we will 
show that the model exhibits saturated ferromagnetic, spin-triplet electron pair condensation 
in the ground state over a certain range of interaction parameters. 

It is worth noting that in the last decade the Majorana edge state formed on a spinlcss 
superconducting wire has attracted much interest both theoretically and experimentally 0EI 
Our model exhibits saturated ferromagnetism where the electrons behave 
as spinless fermions. We show that a similar edge state is formed in the gauge symmetry 
breaking ground state of our model. 

This paper is organized as follows. In Section [21 we give the definition and state the 
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main result. Section [3] is devoted to the proof of the main result. In Section [4l we introduce 
two ground states with broken spin-rotation and gauge symmetries, and show that these 
ground states exhibit similar properties to those of the spinless superconducting wire in the 
topological phase. In Section El we consider mean fields which stabilize the ground states 
introduced in the previous section. In Section [6j we investigate the properties of the ground 
state with the fixed number of electrons. In Section [71 we extend the model to the case of 
the anisotropic spin-spin interaction. Finally, in Section El we provide conclusions. 

2 Definition of the model and the main result 


We consider a one-dimensional array of L sites, which are labeled as 1, 2,..., L. We write A 
for the set of numbers 1, 2,..., L and identify A with the array of L sites. We also write A 
for A\{L}. In this paper L is assumed to be an odd integer with L > 3. This condition is 
adopted only for simplicity, and similar results for even L are obtained with minor changes. 

Let c x ^(c\ a ) be the annihilation (creation) operator of an electron at site x € A and with 
spin a =7, 7- They satisfy the anticommutation relations, 


and 


k 


■X,CT1 Cy,T 


} = k 


,t J 

X ,cn 


} = 0 


( 1 ) 


{ C x,o-y c y,r} ~ fix,yd<j,T ( 2 ) 

for any sites x, y and any a, r =7,7- For each site x, we define the number operators 
= 4,„Cx,„ and », = n x , + and the spin operators = i £ 4X> Cl , with 

(7,T 

/ = 1,2, 3, where pa\ are the elements of the Pauli matrices 


pW = 


0 1 
1 0 


, p {2) = 


0 -i 
i 0 


, P (3) = 


1 0 
0 -1 


(3) 


For each nearest neighbour pair of sites x and x + 1, we define local Hamiltonian H x by 


H x Hi x T Hjj x T Hy x T Hj x T Hx tX 


(4) 
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where 


H t , x 

1 'y ^ ( C x,cr C x+l,cr T C x+l,cr C x,cr ) h-x^x Ax+l^x+l; 

(5) 


cr=t,t 


Hu, x 

U(n x ^n x ^ T , 

(6) 

Hv, x 


(7) 

Hj, x 

7 (n x n x +i c c \ 

~ J V 4 * a+ V ’ 

(8) 

H x , x 

^ ^ (X x Tlx,—cr H - cr) (^ ?cr ^+l,cr ^ * 

(9) 


<T=td 



The term H t x represents electron hopping, and H Ux , H V)X i Hj, x and H x . x represent electron- 
electron interactions, usually referred to as the on-site, the charge-charge, the spin-spin and 
the bond-charge interactions, respectively. In this paper, we assume 0 < 2 1 < V and define 
parameter 6 ranging from 0 to 7t/2 by 

P 2 1 , . 

sm 8 — (10) 

We then consider the Hamiltonian given by 


h = '£h, 

xGA 


( 11 ) 


on A with open boundary conditions. 

Before stating our main result, we have to introduce some more notations. Let us define 
a operators by 


sin 5 


E 

\V =i 


Wy c y,(7 




sin 5 


E 

y =i 


^yCy,cr 


WyCy,a ) if X G A; 

y=X -\-1 

if x — L 


( 12 ) 


where 


J sin(5/2) if x is odd; 
) cos [8/ 2) otherwise. 


( 13 ) 
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By using the a operators, we define pair operators Q T with cr, r =f, | by 




x,y(zA 


(14) 


where F Xiy is given by 


^ sin 8 if y — x = 1; 
-isin<5 ifx = l,y = L; 
F x,y = | sin 5 if x — y = 1; 

| sin 5 if x = L, y — 1; 


(15) 


0 


otherwise. 


It is noted that F x>y = — F yx . 


We denote by $ 0 the state with no electrons on A. The total number of electrons on A 
is denoted by Ay. We assume 0 < N e < L and define the number N p of electron pairs by 


N P = 


N e 


N e ~ 1 


for even N e ; 


for odd N P . 


(16) 


With the values of the parameters given by 

4t 2 


Un — V sin' <5 = 


V ’ 


4 t 2 


J 0 = V(2 - sin 2 5) = 2V- —, 


V 


V 

Xn = —sin 5 cos 5 — t\l 1 


(17) 

(18) 

(19) 


our main result is summarized as follows: 


Proposition 2.1 Suppose that both U > Uq and J > J 0 are satisfied. Then, the ground 
state energy of H with 

X x = (-l) x+1 X 0 (20) 
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and 


/i* = “{l-(-l ) x cos<f} (21) 

is zero for 0 < N e < L. For fixed N e , the ground state is unique apart from the degeneracy 
due to the spin-rotation symmetry, and is given by 

[ (eft) P $o for even N e 

(22) 

( a[ t $o for odd N e 

and its SU( 2) rotations. 

The parameters in the Hamiltonian FI must satisfy several conditions so that $g can 
become the ground state of H. Here we briefly comment on the physical feasibility of 
these conditions. Firstly we need sufficiently large on-site repulsion and nearest neighbour 
ferromagnetic interaction. These are necessary mainly to stabilize the ferromagnetic state. 
(Note that the Hamiltonian H is proved to exhibit metallic ferromagnetism for J > 0 in 
the limit U —> oo [E].) Ferromagnetic materials are expected to satisfy these conditions. 
As for the nearest neighbour charge-charge interaction, it must be attractive. Although 
the charge-charge interaction arising directly from the Coulomb interaction is repulsive, it 
may become an effective attractive interaction such as a phonon-mediated interaction. The 
strength V of the charge-charge interaction also needs to be 2 1 < V. This condition will hold 
in systems with narrow conduction bands. We furthermore need to fine-tune p x and X x . 
Note, however, that the anisotropic spin-spin interaction removes the condition on X x (see 
Proposition 17. II) . Although it will be difficult to fine-tune p x and X x , the ground state $0 or 
a state which has a large overlap with <f>G may be realized in one-dimensional ferromagnetic 
materials with narrow conduction bands if an effective attraction between electrons can be 
generated in the system. See also Section [5] where we treat the superconducting paring field. 
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3 Proof 


Proof of Proposition \2.1l In the following, we assume that the conditions (1201) and (T2Tj) are 
satisfied. We also assume that the electron number N e is fixed. 

Firstly we shall show that the Hamiltonian PI can be expressed as a sum of positive 
semi-definite operators. We define a operators by 




(23) 


for x G A and 


^L,cr ^2^1,cr ~b ^L— l^L,cr* 


We also define b operators by 


(24) 


fox,(T Wx Qc,cr + ^r+lQc+l,cr 


(25) 


for x G A and 


br,cr ^1^1,(7 WLCL ,cr • 


By using the a operators and the b operators we define 


(26) 





(27) 


It is noted that is positive semi-definite. Then, after a lengthy but straightforward 
calculation, one finds that H x is rewritten as 


H x — Hq, x + + Hj\ x + Hw,xi (28) 

where Hu',x and Hj^ x are, respectively, defined by ([6]) and ([8]) with U and J replaced by 
U' = U — Uq and J' = J — J 0 , and H w , x is defined by 

H w,x = ^(4,t4^ + c ]Ui,t4+iJ(^rf^,t + ^+U^+i,t) ( 29 ) 
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with W = Uq/2. For U > Uq and J > Jo, all the terms in the right hand side of (1281) are 
positive semi-definite. This proves that H is the sum of the positive semi-definite operators 
for U > Uo, J > Jo- Therefore, a zero energy state of H, if it exists, is a ground state. 

Let us next show that <f> G is a zero energy state of all the terms in (128|) for any x G A. 
Note that the a operators form a basis for fermion operators on A, since {a\, a , a yj(T } = 
5 XjV for x, y G A by our definition. So we expand b Xj(7 with x G A in terms of a Xi<J as 
b x ,a = Y.y&k{ a l,o, bx,a}a y ,a- It is easy to see that {a\^, b X)fy } = F VtX , which gives us 

bx,cr ^ ^ ^ y,x^y,a• ( 30 ) 

y &A 

From this expression of the b operators we obtain 



Since (fr^) 2 = 0, Q3Tj) implies that commutes with (b\. ^a x ^ + b' x ^a x ^) for x G A. The 
creation operator a\ ^ anticommutes with a x ^, i.e, it also commutes with (b' x ^a x ^ + b\ 
for x G A. Therefore, we have = 0. This together with the fact that there is no 

creation operator with the 4,-spin in <f> G leads to H x & G = 0 for any x G A. This proves that 
<f> G is a zero energy state of H. From the c operator representation of <F G (see Appendix IX|) . 
we find that the ground state is not the null state. 

Finally we shall show the uniqueness of the zero energy state. 

Let M be the eigenvalue of the third component of the total spin. Since the Hamiltonian 
H has the spin-rotation symmetry, it is convenient to decompose the Hilbert space 1~L of 
states into the subspaces 1~Lm each of which has the fixed eigenvalue M. Let <f>M be a lowest- 
energy state in 'Km- Since the representative of <f> G in Km is also a zero energy state of 


H , the lowest energy in T~Lm is guaranteed to be zero. This implies that <3 >m must satisfy 
H x &m = 0 for x G A. In particular, for U > Uq, c x ^c x ^m nrust be zero for any x G A. Now 
we represent ( I>m by using the c operators. As mentioned above, since each site is forbidden 
to be doubly occupied by electrons in it can be expanded in terms of normalized basis 
states in the form 



where A is a subset of A with |A| = N e , a x =t,I, and Yh x &A a ^ = In the product, the c 
operators are ordered in such a way that the site indexes x increase from left to right. 

Let us consider the matrix representation H of the Hamiltonian H with respect to the 
basis states in the form (j32j) . We assume that the basis states are ordered in an arbitrary 
manner and denote by H, :j the matrix element corresponding to i -th and j'-tli basis states. 
Then one easily finds that any non-zero off-diagonal matrix element is —t or —Jj 2, which 
is negative. It is also easy to see that for any i,j there is a sequence A, * 2 , ■ • •, A such that 
Hj^ ... 7 ^ 0. Therefore it follows from the Perron-Frobenius theorem that the 

lowest energy state of H is unique [15], which implies that the lowest energy state of H in 
%m is also unique and is given by the representative of in 'Hm- This completes the proof 
of Proposition 12.11 

Before ending this section, we make a remark on the related exact results of the extended 
Hubbard model. In the above proof we have rewritten the Hamiltonian as a sum of the 
positive semi-definite operators, and then have shown that the ground state attains the 
lowest eigenvalue, zero, of these operators. This strategy was used in Ref. Jb 6 ] to determine 
a parameter range for which the extended Hubbard model has the ferromagnetic ground 
states at half-filling (where the number of electrons is equal to that of sites). Our result 
corresponds to an extension of Ref. [T 6 ] to the case away from half-filling. Note, however, 
that our method for the construction of the exact ground state away from half-filling is quite 
different from that at half-filling. 
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4 Ground states with broken spin-rotation and gauge 
symmetries 

In this section, we assume that the parameters U, J, X x and /j, x satisfy the conditions in 
Proposition 12.11 and hence the ground states of H with the fixed electron number are given 
by (|22l) and its SU(2) rotations. Since the ground states are saturated ferromagnetic, we 
furthermore assume that the third component of the total spin is fixed to N e / 2. In the 
following, since all the electrons are assumed to have the 'f-spin, we omit the spin indexes in 
the fermion operators for notational simplicity. 

The spin-triplet electron pairing ground state of H is regarded as the pairing state of 
spinless fermions. The ground state of our model is thus expected to have some similar 
aspects to that of the Kitaev chain model in which there appears the Majorana edge state 
at the ends of the chain. We will show that it is the case. 

Let us define the zero energy ground states with the broken gauge symmetry 



and 

$ G ,i = V 2 V sin <5 exp (-|e~ ie C f ) $0 = V 2 V sin 5 a[<h G ,o, (34) 

where r/ is a positive parameter and 9 is a phase parameter (note that ("1 = Cft)- The 
state < Lg ) o( < =I?g,i) is a superposition of the zero energy states of H with even(odd) numbers 
of electrons. The states $ 0,0 and 3 >g,i have the different fermionic parities, and, as we shall 
see in the next section, these states are stabilized by superconducting pairing fields. 

As usual, let us define the Majorana fermion operators 

i— —i— + 

7 A,x = e2 Cx + e 2 4, 

7b, x = -i e^c x + ie _1 ^4, 
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(35) 

(36) 





which satisfy = ^/ a>x and {'y a , x ,'Yp l3/ } = 28 at p5 x , y for any a,/3 e {A, B} and x,y e A. By 
using 7 aj i and 7 ^ 1 , with a = A,B } we introduce new edge operators as 


Ti 

- m —t— r {(^2 + Wi)7a,i + 1 (w 2 77 ^ 1 ) 7b, 1 }, 

y/2r] sm 0 

(37) 

r l 

- /7T— T — r {(^2 + V w i)1b,l 1 (w 2 rjWi)i a,l}- 

V zr/ sm 0 

(38) 

(Recall that w± = sin 

(5/2) and w 2 = cos(5/2).) The edge operators Y 1 and Yl 

are rewritten 

as 




T, = */— ^—-(w 2 e 1 ^ci + r]w 1 e~ 1 ^c\') 

(39) 


V Tj sm 0 V / 


r L = 1 a / . J. ( w 2 e 1 ^c L + r]w 1 e 

V 7/sino V / 

(40) 

with the c operators. 

Then, we find that 



ri$ G;0 = — irL<F Gj o = e^ l2< F Gi i, 

(41) 


ri$ G ,i = ir L $ Gjl = e l2 $ G0 . 

(42) 


Furthermore, from the above relations, we obtain 


— irir L <f>G,o = <h G ,o, (43) 

—irir L $ Gil = —$c,i- (44) 


The relations (TUI) and (1421) are obtained as follows. For 16 Awe have from (l3Tf that 


i*a x (-|e- ie C f ) n $0 = ve-'Mn (~f e" ie C f )” ' $0, 


which yields 


$o n = r?e ' 2 bl <3> 


e 2 a x <J > G)0 = V e 


G,0- 


(45) 


(46) 


Here we used b\. ((4) 2 $0 = 0 [ITT] . By (146p . we also have 

e l2 a x $ G ,i = e l2 a x yy/2r]am8 a\^ G ,o) = ^e - ' 2 ^.*f> G ,i + 5 XtL y/2r)B\\\5e 1 *$ G ja- (47) 
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By representing (146|) and (|47)) with the c operators and setting x = L, one finds 


{ w 2 e^ (ci + c L ) + r]Wie *2 - 4) } $g ,i = Si ,i y/2rismS e 1 ^ $ G ,o (48) 

with l — 0 , 1 . On the other hand, (ITS]) and (147T) combined with 

b\. = —w\c\ — wic\ + 2 sin 5a) L (49) 

xgA 

y ^g x = VJ 2 CI-W 2 C L , (50) 

x£A 

which follow from the definition, yield 

{ w 2 e^ (ci - Cl) + rjw (4 + 4) | = S ifi \/2ri sin5e _ 1 ^<f) Gi i (51) 


with / = 0,1. From (1481) and (1511) we obtain (jUJ) , and (J42]). 

It is noted that in the case 7 = w 2 /w± = 1/tan(<5/2) we have F] = 7 and Yl = 7 b,l 
which are the Majorana fermion operators. In this case, we can reconstruct the edge fermion 
operator by combining 7 a,i and 7 b,l as 


d e d g e = \ e l ° 2 (7a, 1 + 17b, l) ■ 

The fermion operator d e d g e satisfies {d e d g e, d edge } = (4d ge > 4d ge l = 0 and (4d g e> ^edge} 
From (EXT]) and (EXT]) we also have 


(52) 
= 1 . 


4dge$G,0 - $G,1, (53) 

d e d g e$G,l = ^*G,0- (54) 

The above relations yield n edge $ Gjl = $ Gjl and ri edge $ Gj o = 0 with n edge = d\ Age d e d g e, which 
imply that the Majorana edge state is formed at the ends of the chain. 


5 Mean field Hamiltonian 

In this section we consider external fields (or mean fields) which remove the ground state 
degeneracy and select $ G o and $ Gj i as the two ground states. 


12 




It is well known that the external magnetic field can remove the degeneracy due to the 
spin-rotation symmetry. So we assume that the system is in a magnetic field, and fix the 
third component of the total spin to N e /2. (As in the previous section, the spin indexes are 
omitted in this and the next sections under this assumption. ) 

In order to remove the degeneracy due to the electron pair condensation, we shall consider 
the Hamiltonian which does not conserve the electron number. More precisely, we will 
introduce Hamiltonian H' of spinless fermions with superconducting pairing field, and show 
that the ground states of H + H' are given by $ 0,0 and $g,i- 
Let us define 


H' 

Hi 


/ 


2> 

xGA 

A (e -> 

n 


(e '*4 -rie l °b x )(a + (1 - a)a x al)(e'°a x -r]n 


(55) 

(56) 


where a and |A| are non-negative parameters. As we will see below, A = |A|e l£l corresponds 
to the superconducting pairing field. Since {aj., a x } = 1 for x G A, we have 


Hi = 


V 


(e ‘=4 - - (1 - a) 4 aJ(e' 2 o x - 7)e ’=), 


(57) 


and hence H' x is a positive semi-definite operator for a > 0. From (|46|) and (147|) we find that 
<Fg,o an d $gp are zero energy states of H' x for x G A. Therefore $g,o and $g,i are ground 
states of H + Hi It is easy to see that there is no other ground state. The Hamiltonian H' 
removes the ground state degeneracy of H and stabilizes the states <3>g,o and <I>Gp- 
After some lengthy but straightforward calculations, H' is rewritten as 


H' 


with 


S 'y ^ { C x C x+l T C x+l C x) ^ T ^x+l C x+l C x+l) 

ieGA x£A 

- v '^2 c lcxci +1 c x+ i + ^2 (Ac x c x+1 + A*4 +1 4) + rj\A\(L - 1) 

xGA xGA 


s = 


|A| 

2r] 


(1 + arj 2 ) sin <5, 


(58) 


(59) 
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(60) 

( 61 ) 


14 = {(1 + aij 2 - 2i; 2 ) - (-1)*(1 + oij 2 ) cosi} , 

V' = (a-l)i;|A|. 

From the above representation of H', one immediately realizes that A corresponds to the 
superconducting pairing held, which may be induced from a nearby superconductor. This 
held term essentially removes the degeneracy. It is noted that, in the case where 5 = ir/2, 
rj — 1 and a = 1, H' is reduced to the Hamiltonian of the Kitaev chain of the spinless 
fermions in the topological phase. Thus our model can be also regarded as an extension of 
the Kitaev chain to the spinful system with the electron-electron interactions. 

6 Electron number conserving case 

In the previous two sections we considered the case where the number of electrons is not 
conserved. From the expressions (fool) and (154)) of the symmetry breaking ground states, one 
hnds that the edge state is closely related to the zero energy mode corresponding to a) L . 
Indeed, we have shown that the occupation of a) L by an electron is reflected as an eigenvalue 
of the number operator n e d ge of the edge fermion operator. 

For the fixed electron number, the ground state <f>G can not be the eigenstate of n e d ge , 
since we have 

^ed g e = ^ (1 + i7A,i7B,n) (62) 

with 

i7A,i7B,L = e'%c L + e~ l6 c[c\ + c\c L + cjci. (63) 

Instead, we can expect that there is a difference between the expectation values of n e d ge for 
<Fg with N e even and odd. 

Let (• • -) 0 and (• • -)i be the expectation values (<I>g> • ■ • $g)/^g) f° r ^g with N e even 
and odd, respectively. We will estimate (n e d ge )o an d (^ e d ge )i- Clearly, we have {ciCl)i = 


14 


(c^ L c\)i = 0 and (c[cl)i = (cj^ci)z with / = 0,1. Let us consider ( c\cl)i ■ By using the c 
operator representation of $g (see Appendix [Aj, we obtain 


(c\c L )i = (-1) 


1+1 sin 2 I - 


XMcA;|A|=AT e -l XWi L ^ A] W A 
2 / ^AcA;|A|=iV e Wa 


(64) 


where Wa = IIxeA w x> an d x[E] takes the value 1 if E is true and 0 otherwise. Since we 
have 


X] W A < 

AcA;|A|=V e 


L(L-l) 

N.(L - i V e ) 


cos 2 


]T X [1 ,L£A]W A 

AcA;|A|=AT e -l 


(see Appendix iBlh \( c [ cl ) i \ is bounded from below as 

\(c\c L )i\ > tan 2 ^0 p(l-p) 
with p = N e /L. Therefore, we obtain 

0 < i - tan 2 (( ) fit I - p) 
i > i + tan 2 ((j fill - p). 




(65) 


( 66 ) 


(67) 

( 68 ) 


The inequalities obtained above show that the ground state expectation value of the 
occupation number n e d ge corresponding to the edge fermion reconstructed by the Majorana 
fermions depends on the fermionic parity, regardless of the chain length L. Let N e be even. 
Since we have p = N e /L & (N e + 1)/L & (N e — 1)/L for sufficiently large L, these inequalities 
imply that the expectation value of n e d ge for N e decreases by at least 2 tan 2 (<5/2)p(l — p) 
compared with that for A’ e — 1, while the expectation value of n e d ge for N e + 1 increases by 
at least 2tan 2 (5/2)p(l — p) compared with that for N e . This behavior in the edge fermion 
number indicates the formation of an edge state in the electron number conserving setting. 
In the following, we propose a concrete example of a system having the two-fold degenerate 
ground states each of which is characterized by a zero energy mode related to the Majorana 
edge state. Very recently, a similar model has been investigated in Refs. [jT8] and [ T9] . 
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Firstly we prepare a copy of H. The operators in the copied system are denoted by the 
underline as c x . We then consider the Hamiltonian H + H_ + H e on the two chains, where 






(69) 


with e > 0 is an interchain interaction. The number of electrons on the whole system 
is fixed to N e . We suppose that the values of the parameters in H and H_ are taken so 
that each Hamiltonian is positive semi-definite and has the zero energy ground states (see 
Proposition 12.ip . Under the assumption that the system is in a magnetic field, one finds 
that the two states 


*o,o = “l (C' + C , )' Vp 4 > ci, *0,1 = “1 (< t + C , ) Wp * 


(70) 


for odd N e , and 


*o,o = (C* + C t ) Wp *o, *y =aULU + C ')' 


4>f 


(71) 


for even N e are the only ground states of this system. In fact, //. H_ and //, are positive 
semi-definite, and $ G 0 and <I> G j are the only zero energy states for these Hamiltonians. 

It is expected that similar inequalities corresponding to (16711 and (l68p hold for the ground 
state expectation values of the number operators n e d ge and n ed ge of the edge fermions on the 
chains, although explicit analytical expressions are difficult to obtain. 

We end this section with the remark that the fermion operator defined by a n = XLga( — 1) x+1 Ox 
plays an interesting role in manipulating the zero energy mode in the condensate. More pre¬ 
cisely, a,satisfies since XLeA( — l) x+1 ^ = 0 and {a^ L , a n } = {a^, a^} = 1. There¬ 

fore, we have the relations {yJ2rj sin 5) _1 a 7r <l>G,i — ^g.O; ^L a ^G,i — ^g.i and ^L a n&G,o — 0 
for the symmetry breaking ground states. Similar relations are also found for the elec¬ 
tron number conserving system. Namely, we have a^a T $' G0 = $ G0 , a} L a n & G 1 = <P G x and 
o = 4 a > t$g,i = 0 for odd N e , and a{a n $ GA = G1 = <P G1 and a\a n <$>' Gfi = 

a L<2 7 r < 3?G o = 0 for even N e . 
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7 Spin-Spin Interaction with Ising-like Anisotropy 


In this section, we treat the case of the spin-spin interaction with an Ising-like anisotropy. 
Let us define 

(S* ■ s x+ 1 )„ = + P(SVS™ 1 + Sf sg,) (72) 

where f3 is a non-negative parameter and denote by Hjp jX the Hamiltonian obtained by 
replacing S x ■ S x+ 1 with (S^ • in H Jx . Then we consider the Hamiltonian 

h p = ( 73 ) 

xGA 

Hp,x = H tx + H Ux + H Vx + (74) 

Note that the bond-charge interaction H Xx is omitted in Hp )X . For Hp, we have the following 
result: 


Proposition 7.1 Suppose that both U > Uo + 2A" 0 and J > Jo + 4X 0 are satisfied. We 

furthermore suppose that p, x is given by / OH]) . Then, the ground state energy of Hp with 

4X 0 


T<p<i 


J 


is zero. For fixed N e , the ground state is two-fold degenerate and is given by 


_ J (CUF" to for even N, 

/oraMJV. 


■I',: = 


(75) 


(76) 


with a =t, 4" 


The outline of the proof is as follows. As in the isotropic spin-spin interaction case, we 
rewrite Hp x as 

Hp, x — Ho,x + Hu",x + Hj",p\ x + H x 0 ,x + H Wx , (77) 

where H Xo , x is given by 

Hx o,® = XoY^ {4,«r + (-l)*4+i,«x} n x - a {c x , a + (-1 ) x c x+ha } 

o-=t,4 

+Xo { 4 , a - (-l) X 4+i I<r } W*+l,-«T {Cx,a - (-1)*^+!^} , (78) 

c=t,4 
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Hu", x is obtained by replacing U with U" = U — Uq — 2X 0 in Hu tX , and Hjn^ x is obtained by 
replacing J and (5 with J" = J — J 0 — AX 0 and /3 1 = (J f3 — Jq) / (J — Jo — ^X 0 ), respectively, in 
H JAx . When U" > 0, J" > 0 and 0 < /3' < 1, all the terms in (1771) are positive semi-definite 
and in (1761) is their zero energy state. The fact that there is no other zero energy state 
follows from the application of the Perron-Frobenius theorem. 

In the case of the isotropic spin-spin interaction, the bond-charge interaction whose 
strength parameter is fixed must be included in the Hamiltonian to obtain the exact ground 
states. On the other hand, the Hamiltonian with the anisotropic spin-spin interaction has 
the exact ground states even if the bond-charge interaction is absent. Although the on-site 
potentials still have to be adjusted to certain values, the model with the anisotropic spin- 
spin interaction exhibits the spin-triplet electron pair condensation over the wide range of 
parameters. 

8 Conclusion 

We have introduced the one-dimensional extended Hubbard model whose ground state si¬ 
multaneously exhibits saturated ferromagnetism and spin-triplet electron pair condensation 
under certain conditions. Recently, the extended Hubbard chain with charge-charge and 
spin-spin interactions at low filling was studied by means of mean field and numerical meth¬ 
ods in Ref. [SO] . The results showed that the ground state is in the spin-triplet pairing phase 
for strong ferromagnetic coupling, even if there are no fine-tuned bond-charge interactions 
and on-site potentials which are necessary to get our exact results. These results together 
with ours indicate that the model exhibits spin-triplet pairing over a wide range of param¬ 
eters. We have constructed two ground states in which both of the spin-rotation symmetry 
and the gauge symmetry are broken. It has been shown that these ground states are trans¬ 
ferred from one to the other by applying the edge operators. The edge operators become 
the Majorana fermions in a certain case, and, in this sense, the Majorana state is formed on 
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the edges of a chain in our model. We have introduced the mean field Hamiltonian with the 
pairing held which stabilizes the gauge symmetry breaking ground states. Here we remark 
that the spin-triplet pair condensation found in the ground state of H is unstable against 
the thermal fluctuation since H is constituted of short-range interactions and is defined on 
a chain. However we can expect that the spin-triplet pair condensate survives at non-zero 
temperatures in the strong pairing held. The mean held Hamiltonian has been shown to 
be regarded as the Kitaev chain with the nearest neighbour charge-charge interaction. It 
is noted that a similar spinless fermion model has been studied by Katsura, Schuricht, and 
Takahashi recently [21]. Our extended Hubbard model together with the mean held is an 
extension of the Kitaev chain to the spinful electron model. We have also estimated the 
expectation values of the edge fermion number operator for the ground states with hxed 
even and odd numbers of electrons, and found that there is the difference between them. 
Furthermore, we have proposed the model on the two chains in the electron number con¬ 
serving setting and have shown that the model has the two-fold degenerate ground states 
which are characterized by the zero modes on the chains. 

To conclude, it is interesting to note that Nadj-Perge et al. reported the observation of 
Majorana fermions in a chain of Fe atoms, which intrinsically have ferromagnetic nature, 
on a superconducting Pb substrate MM- It is also noted that the recent developments 
in the field of cold atoms open a route to the experimental realization of one-dimensional 
interacting fermion systems [22] . We hope that our results stimulate these fields. 
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Appendix A c operator representation of <£>g 


Here we express <J>g in terms of the c operators. For the notational simplicity, we omit the 
spin indexes of the fermion operators. 

From (TT5T) and (ITlj) . one obtains 


d = — sin 


in5 

\xGA 


® x a x+i + a\a} L | . 


Substituting (fl2l) into the right hand side of (1A.1I) . we have 


c f = 


4 


sin 5 


54 0 ^ 44 )- 


(A.l) 


(A.2) 


x,y€A,x<y 

Then, taking into account the sign factor arising from the exchange of fermion operators, we 
have 


and 


(4)^0 = 


4 (4)^0 = 


4 \*p 


sin 5 


W) 


54 n^ 4 ) 


i AcA;|n|=2W p xeA 


4 \*p 


sin 5 V sin 5 


W) 


54 n^ 4 ) 


(A.3) 


(A.4) 


, J 4cA;|A|=2A r p+l x&A 

where N p \ = N P (N P — 1) • • • 2 • 1, and | A\ denotes the number of elements in a set A. 


Appendix B Proof of the inequality ( 1651 ) 


Let us prove the inequality (J65J). Firstly we rewrite the left hand side as 


E w * = wY, w * E x\x£A]w a . 


N e A 

AcA;|A|=W e x&A AcA;|A|=We-l 

Then, by using Wl = sin(<5/2) < w x < cos(<5/2) (recall that 0 < <5 < vr/2), we obtain 

54 Wa - f cos2 (T)54 54 x[x^a}w a 

AcA-,\A\=N e e ' ' x&A AcA;|A|=We-l 


(B.l) 


5 K cos 


E x[L £ A] W A . 


(B.2) 


AcA;|A|=V e -l 
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Here, note that x[L ^ A] — %[1 ,L ^ A] + x[l G A, L ^ A], Since the sum related to 
x[l G A, L £ A] is bounded as 


E x[l e A,L £ A]W a 

AcA;\A\=N e -l 

A E x[1,L£A]^a 


= sin" 


AcA^yl^Afe—2 


5 

- sin | - 




w; 


AcA;|A|=A r e -2 


ce A 


w x(L ~ N e ) 


X[x £ A U {1, L}} 


< 


L — N e 

N e - 1 
L-N P 


E E x[l, L ^ A]y[rr A U {1, L}] 


AcA;|A|=A r e -2 zeA 


E 


(B-3) 


AcA;|A|=N e -l 


we conclude that 


E ^ - at ( 1 + 


AcA;|A|=JV e 


iVe- 1 
L — N P 


cos - 


E x[l ,L$A]W a 


AcA;|A|=N e -l 


L- 1 


TVp \ L — N P 


cos - 




(B.4) 


AcA^A^Ve-l 
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